Density-functional formulations for quantum chains 
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We show that a lattice formulation of density-functional theory (DFT), guided by renormalization- 
group concepts, can be used to obtain numerical predictions of energy gaps, spin-density profiles, 
critical exponents, sound velocities, surface energies and conformal anomalies of spatially inhomo- 
geneous quantum spin chains. To this end we (i) cast the formalism of DFT in the notation of 
quantum-spin chains, to make the powerful methods and concepts developed in ab initio DFT avail- 
able to workers in this field; (ii) explore to what extent simple local approximations in the spirit 
of the local-density approximation (LDA), can be used to predict critical exponents and conformal 
anomalies of quantum spin models; (iii) propose and explore various nonlocal approximations, de- 
pending on the size of the system, or on its average density in addition to the local density. These 
nonlocal functionals turn out to be superior to LDA functionals. 
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I. INTRODUCTION 

In this work we show that a lattice formula- 
tion of density- functional theory (DFT), guided by 
renormalization-group concepts, can be used to obtain 
numerical predictions of energy gaps, spin-density pro- 
files, critical exponents and conformal anomalies of spa- 
tially inhomogeneous quantum spin chains. 

Density functional theory emerged over the years as an 
efficient and powerfuU tool for electronic-structure theory 
in condensed matter physics, material science and quan- 
tum chemistry. In spite of these successful applications 
of DFT to continuum systems, ^'^ DFT is not much used 
for interacting discrete systems, i.e., systems where the 
operators are defined on a lattice. Previous applications 
of DFT for discrete systems were to the one-dimensional 
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or 



Hubbard model describing electrons in solids 
atoms in optical latticesj^i^^ and to the Heisenberg spin- 
i model in d-dimensionS )^^'^^i^^'^'* employing local ap- 
proximations similar to the local-density approximation 
(LDA) of ab initio electronic structure theory. The diffi- 
culty in treating such systems resides in the construction 
of a good functional. In Ref. i, an LDA was constructed 
for the Hubbard model by exploiting the exact Bethe 
ansatz solution for the ground state of the infinite, spa- 
tially homogeneous, system. In Ref. [ll|, LDA-like func- 
tionals for the Heisenberg model were obtained from an 
approximate spin-wave theory for the infinite homoge- 
neous system, improved by adding numerical corrections 
from density-matrix renormalization-group calculations. 

The purpose of this paper is to introduce a possi- 
ble DFT formulation appropriate for general quantum 
chains, in particular near their quantum critical points. 
The construction of approximate density functionals ex- 
ploits two seminal concepts from the theory of critical 
phenomena: the renormalization group and the confor- 
mal invariance of critical systems. From the renormaliza- 



tion group, it is known that the critical systems can be 
classified in universality classes. In each of these classes 
the low-energy physics is described by distinct homoge- 
neous systems (fixed points) with correlation functions 
defined by appropriate critical exponents. Moreover, the 
underlying field theory describing most of these univer- 
sality classes is conformally invariant. This symmetry 
classifies these classes by the conformal anomaly c, or 
the central charge of the conformal algebra, and the crit- 
ical exponents are given in terms of the highest weights 
representations of this algebra. These facts lead us to 
expect that a single density functional could describe the 
whole family of critical models of a given universality 
class of critical behavior. In principle, the construction 
of such functionals can proceed by exploiting the exact 
results provided along the years by the Bethe ansatz^'' in 
its several formulations. 

It is believed that for each universality class of critical 
behavior there is at least one exactly integrable quan- 
tum chain. The exact solution of that quantum chain 
underlies the construction of the functional for several 
models belonging to the universality class. In order to 
test this general idea we restrict ourselves to one of the 
most important universality classes of critical behavior, 
namely, the c — 1 conformal invariant models, with con- 
tinuously varying critical exponents. Models in this uni- 
versality class are the Gaussian modelfii the anisotropic 
antiferromagnetic Heisenberg model with half-odd in- 
teger spin,i^ii2. the ferromagnetic version of the exact 
integrable Taktajan-Babujian models^^ the anisotropic 
biquadratic spin 1 chain, the Ashkin- Teller quantum 
chain, etc. 

As a representative of this class of models we are go- 
ing to consider the exactly integrable spin-i anisotropic 
Heisenberg model or the XXZ quantum chain. This 
model is considered a paradigm of integrability, and its 
critical exponents are known exactly. 
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Due to the conformal invariance, the leading finite-size 
corrections of the eigenenergies of the quantum chains 
are ruled by the conformal anomaly and the critical ex- 
ponents of the critical chains i^^'^'^ The exact knowledge 
of the ground state energy in the bulk limit and in the 
finite geometry will allow us to construct good LDA ap- 
proximations for the functionals. Some of these LDAs 
will be good enough to reproduce the exact values of the 
conformal anomaly and critical exponents. Such LDA 
approximations, although applied in this paper only to 
the XXZ chain, reveal interesting general features, that 
certainly will be helpful for other discrete quantum sys- 
tems. Our analysis reveals the essential ingredients and 
the degree of accuracy we should demand from a den- 
sity functional to obtain from it reliable results for the 
conformal anomaly and critical exponents. 

Once such a functional is available, numerical analy- 
sis of complex Hamiltonians is greatly simplified. As an 
example, in the case we are mainly concerned with be- 
low, the XXZ chain, full exact diagonalization becomes 
hard for more than « 12 sites, Lanczos methods allow 
one to reach about three times as much, whereas density- 
functional techniques can still be applied for systems with 
thousands of sites, even if translational invariance is bro- 
ken. 

The layout of this papers is as follows. In Sec. [Til we 
give an introduction to the general framework of DFT 
for spin chains (Sec. Ill A[) and to local approximations 
that make this framework useful in practice (Sec. IIIB[) . 
In Sec, mil we then describe the representative spin chain 
we are mainly concerned with in this work: the XXZ 
model. Section IIVI is devoted to the construction of local 
approximations for the XXZ chain, and an exploration of 
their performance for systems with periodic fSec. IIV"S]) . 
twisted (Sec. IIVBI) and open (Sec. IIVCP boundary con- 
ditions. In Sec. |V] we then go beyond local approxi- 
mations, and construct nonlocal functionals (Sec. IV A[) . 
which are applied to systems with open boundary condi- 
tions (Sec. IV Bl) and in the presence of spatially varying 
magnetic fields (Sec. IV C|) . Section IVTl contains our con- 
clusions. 



II. DENSITY-FUNCTIONAL APPROACH TO 
QUANTUM SPIN CHAINS 

A. General aspects of DFT for quantum spin 
chains 

Density-functional theory is a formally exact way to 
cast the many-body problem in terms of densities instead 
of wave functions. For a large class of Hamiltonians — 
characterized by a stable ground state and containing a 
density-like intensive variable coupled to an external field 
— the fundamental Hohenberg-Kohn theorem shows that 
the density variable determines the ground state wave 
function, which implies that all ground state observables 
are functionals of the densityjii^ 



This general theorem does not tell one how to obtain 
these functionals, how to calculate the density, what the 
physical nature of the external field is — and not even 
how to select the density variable in the first place. All of 
these questions must be answered before applying DFT 
to a specific class of problems. Here we attempt to do 
this for a large family of spin Hamiltonians, by combining 
concepts from DFT with renormalization-group ideas. 

In practice, the density variable is typically chosen to 
represent the charge or spin distribution in a continuous 
space or on a lattice. The ground state energy becomes 
a functional of this distribution, which is minimized to 
provide the ground state density. Of course, the energy 
functional is not known exactly for nontrivial many-body 
systems. Successful approximation schemes are based on 
the local-density concept, which takes as an input the 
ground state energy of a spatially homogeneous system, 
in which the density distribution is uniform. The corre- 
sponding energy density is then evaluated point by point 
at the actual densities of the inhomogcneous system un- 
der study. 

While this means that DFT in the local-density ap- 
proximation (LDA) requires as an input a good solution 
to the spatial homogeneous system (which in itself may 
be hard to obtain), it also implies that once that solution 
is available, the LDA can be used to obtain approxima- 
tions to the ground state energy and density distributions 
of inhomogcneous systems. In ab initio applications of 
DFT, the homogeneous system is the electron liquid, the 
external potential is produced by the nuclei, and the re- 
sulting inhomogcneous systems are atoms, molecules and 
solids, 

For model Hamiltonians defined on lattices, the ho- 
mogeneous system is one in which all sites are equiva- 
lent, and inhomogcneous systems can have boundaries, 
impurities, confining potentials, or other translation- 
symmetry breaking terms. Since model Hamiltonians are 
typically much better controlled than the ab initio one, 
much more information is available for constructing the 
LDA and for testing its predictions. This additional in- 
formation provides a two-pronged opportunity: (i) Learn 
about inhomogcneous model Hamiltonians, which have 
too few symmetries to be integrable and require too large 
matrices to be numerically diagonalizable, by means of 
LDA calculations, (ii) Learn about DFT, the LDA, and 
its improvements, by studying them in the context of 
model Hamiltonians. Here we explore both aspects for 
spin Hamiltonians, and inquire, in particular, if and how 
LDA-type functionals can be used to predict the critical 
exponents and conformal anomalies associated with the 
quantum critical points of these Hamiltonians. 



B. Local density functionals for quantum spin 
chains 

The family of quantum chains for which we are go- 
ing to seek a DFT formulation, describes quantum spin 
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operators attached to L sites (i = 1, . . . , L) of a one di- 
mensional lattice, with Hamiltonian 

H = Ho + V'''. (1) 

The operator Hq may contain, e.g., the kinetic energy op- 
erator and the static interaction energy among the spins 
on distinct sites, while F*"^* is the on-site potential. 

Here, we consider Hamiltonians of the form ([T]) that 
commute with the global charge N = X^^^^ni, whose 
eigenvalues take a discrete set of values {N = 0, 1, 2, . . .). 
Moreover the on-site potential in (jT]) is a function of the 
local density operators n^, having the general form 



i=l 



(2) 



where (i = 1,...,L) are arbitrary site-dependent 
couplings that fix the profile of the spacial inhomo- 
geneities of the quantum chain. Examples of models in 
this family are the spin 5* Heisenberg models, where N is 
related to the S'^-magnetization, and the Hubbard?^ and 
t-J models, where N is related to the number of fermions. 

For this class of models we can separate the Hilbert 
space associated to (HJ into block disjoint sectors la- 
belled by the eigenvalues TV = 0, 1, 2, . . .. The variational 
principle, restricted to a given sector of total charge 
N = X^iLi '^i; gi'^s US the lowest eigenenergy on the sec- 
tor 

Eo = min|^> < >=< M/f^ |i/ 1^^^' >, 

(3) 

where I^'q^^ > is the eigenfunction corresponding to this 
lowest eigenenergy. 

The minimization procedure of Q can be done conve- 
niently in two steps by following the constrained-search 
approach of Levy^^ and Lieb.^^ In the first step we con- 
sider, for a given density distribution rii (i = 1, . . . , L), 
restricted to = ensemble of states {l^* > 

}ni,...,nL such that < >= m {i = 1,...,L). On 

this ensemble we minimize < >, i.e.. 



min|*>„^ „^ <^\H\^>=FI[n]+Y,v^'n^, (4) 



where the special form of ^ was used, and 

Fi[n\ = niin|*>^^ < ^'|i?o|* > (5) 

is in general a functional of the density n. On a lattice, 
Fi[n\ reduces to a function of the L independent vari- 
ables ni, . . . but for consistency with the ab initio 
literature we continue to employ the expression 'func- 
tional'. For notational convenience, we below frequently 
employ n in the argument of a functional, i.e., between 
square brackets, to denote the entire discrete distribu- 
tion of values ni, . . . ,ni,. When used on its own, or as 



an argument of a simple function of one variable, n rep- 
resents just one value, which can be site dependent (ui) 
or constant {n ~ N/L). 

In the second step we minimize Q over all densities 
satisfying J2i=i '^i = i-^-, 



i=l 



(6) 



where /x is a Lagrange multiplier and plays the role of a 
chemical potential. We have then reduced the variational 
problem ([3]) to ([5]) and ©. Unless we have a good ap- 
proximation for this procedure is, of course, purely 
formal. 

In order to proceed, we now suppose that we can de- 
fine a simple Hamiltonian Hq acting on the same Hilbert 
space as ((l} and having the following three properties, (i) 
It has the same global symmetries as Hq in ([T]). (ii) We 
are able to calculate its ground state even in the presence 
of external inhomogeneous on-site potentials of the type 
given in ([2]), where the Hamiltonian is given by 



i=l 



(7) 



with u| (i = 1, . . . , L) arbitrary, (iii) There exists a spe- 
cial choice of ?;| (i = 1, . . . , L) such that the ground state 
of d?]) has the same density {m} as the ground state of 
the nontrivial interacting Hamiltonian ([T]) . The auxiliary 
Hamiltonian Hg is the Kohn-Sham Hamiltonian of DFT. 
It is not intended to approximate the complex problem 
posed by the Hamiltonian H, but to obtain the ground 
state densities and energies of H}'^ 

Property (i) is not rigorously necessary, but it simpli- 
fies our analysis without too much restricting its gen- 
erality. It allows us to consider Hg on the sectors la- 
belled by the eigenvalues N = 0,1,2,... of the com- 
muting operator N. Property (ii) is needed to make 
the formulation practical, as it ensures that we can cal- 
culate the ground state energy and density of the sim- 
pler Hamiltonian ([7]). Property (iii) is equivalent to a 
v-representability assumption,^ and is essential for the 
procedure to work. 

Given these three properties, we can use the variational 
procedure ([2])-® to obtain the ground state energy from 
the minimization 



i=l 



where, as in ([5]), 



Fg [n] — mini 



l*>..i, 



(8) 



(9) 



and /i*' is a Lagrange multiplier. Comparison of ([8]) and 
implies that the density profile (rii,...,ni) corre- 
sponding to the ground state of the complicated Hamil- 
tonian (III) is the same as that of the ground state of the 
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simple Hamiltonian 7?^, provide we choose 

L L 



vtn, = ^r*"* + W[n\ + {^Ji- n')N, (10) 



i=l 

where 

This implies 



W[n]^Fi[n]-F,[r 



dW[r 
dm 



(12) 



For a given site i, the potential w| depends on the entire 
distribution of local densities m, . . . , n^. We have then 
to solve ([7]) selfconsistently for the ground state energy 
Es and the corresponding densities The 
difference /i — /i** in pH)) and is just a harmless con- 
stant in the above sclfconsistcnt procedure, since we are 
always in a sector with a fixed value of N, and hence we 
may choose fi — fi'^ . As a consequence of (|TOl) and ([S])-©, 
the ground state energy of the complicated Hamiltonian 
© is 



A dW[n] 
1^— 



i=l 



(13) 



which can be calculated once the density is known. 

The success of this DFT formulation will depend on 
our ability to produce good approximations for H^[rt]. 
In ab initio DFT, the functional Fs is just the kinetic 
energy of the noninteracting Hamiltonian Hs, and the 
fimctional Fj[n\ is decomposed as^ 



Fi[n] = T,[n] + Enin] + E,,[n], 



(14) 



where Eh is the Hartree energy (the mean-field approx- 
imation to the Coulomb interaction energy) and Exc the 
exchange-correlation energy. E^c is unknown, and must 
be approximated. Eh is known explicitly, and Ts, al- 
though not known explicitly as a density functional, is 
easily obtained from solving the Hamiltonian Hs- For 
general model Hamiltonians, the various terms in Fj need 
not have the same physical interpretation. We therefore 
simply write Fj[n\ — Fs[n\ +W[n\, as above, and develop 
approximations for M^[7i]. 

If Fj contains pieces that are known exactly, it is ad- 
vantageous to treat these separately, and approximate 
only the remainder (e.g., E^c in the ab initio case). A 
typical case is the mean-field approximation to the inter- 
action energy, leading to a Hartree-like contribution to 
W. Below we develop local approximations to but 
in the models we are concerned with the Hartree term it- 
self is local, so that nothing is gained by giving it special 
treatment at this stage. 

In order to turn the formulation into a practical tool, 
we require approximate expressions for W[n\ that keep 
the essential ingredients of the exact unknown functional. 



A general prescription for systematically producing such 
approximations is to split the functional W[n] into L 
identical parts, according to 



) = E 



Wi 



(11) where 



W{ni, ...,nL) 
L 



i = l. 



(15) 



(16) 



We now expand, for a given inhomogcncous distribution 
rii, . . . , n^, each of the L terms Wi around a distinct ho- 
mogeneous distribution. 



3=1 1=1 

The i'th term is thus expanded around a homogeneous 
distribution in which all sites have the density Ui of the 
i'th site in the inhomogcncous distribution. Here aj [ui] 
and /3j,/[n] depend on the local densities. 



1 dW{ni,...,nL) 
L drii 

1 d^W 



2L duidni 



(18) 



The first term in this expansion provides a local density 
approximation (LDA) to the functionals. The second and 
higher-order terms yield non-local contributions similar 
to those obtained in gradient expansions of the density 
functional for Coulomb- interacting electronic systems. 

Most of our analysis in the next sections is based on 
the LDA, i.e., the first term in HT]) 



W 



LDA 



(ni 



i=l 



L 



(19) 



The exact functional W evaluated at an infinite and ho- 
mogeneous density distribution ni = ••• = = n^, 
reduces to the functional of the infinite homogeneous sys- 
tem, depending only on the variable n^, whose value per 
site we denote as w°°(n). The LDA then assumes the 
more familiar form 



W 



LDA.( 



(20) 



To compare this with the LDA in use in electronic- 
structure theory, recall that there the Hartree term is 
nonlocal, and the local approximation is only applied to 
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the difference E^c = Fi — Fs~Eh, whereas here we apply 
it directly ioW = Fi-Fs. 

In this paper we explore the possibility to obtain 
the critical exponents and conformal anomaly of criti- 
cal Hamiltonians, by means of simple functionals con- 
structed following the above prescription. As a conse- 
quence of conformal invarianccr^ these quantities are 
obtained from the amplitudes of the 0(1/L^) finite-size 
corrections of the energies per site of the quantum chains 
of size L. A finite-size quantum chain can still be ho- 
mogeneous, if all sites have the same spin, and periodic 
boundary conditions are applied. Although the LDA be- 
comes exact in the infinite homogeneous system, in a 
finite homogeneous system it will have an additional er- 
ror. If such error is of 0(1/L^), we are going to obtain 
wrong critical exponents. 

An improvement is obtained by using the ground state 
energy per site of a translationally invariant Hamiltonian 
with L sites, instead of that of the infinite system. In this 
case the functional depends on the size L of the system, 
and can be written as 

L 

W^^^-^{n,, . . . , ni) = ^ ^«(^)(n)|„^„., (21) 

i=l 

where e'^^^ni) is the ground state energy per site of the 
Hamiltonian iJp- 

Note that this LDA- like expression is not, strictly 
speaking, a local functional, as it depends on the size 
of the system (i.e., the number of elements in the set 
{ni . . .n^}), which is a highly nonlocal function of rii. 
The resulting effective magnetic field hi at site i does not 
depend only on the density at that site (as it does in local 
approximations), but on the number of sites. In applica- 
tions reported below, this nonlocality endows the func- 
tional with superior properties, compared to the usual 
LDA. In spite of this nonlocality, we call this latter func- 
tional the finite- size LDA, because of its conceptual sim- 
ilarity with the usual (infinite-size) LDA, deriving from 
the thermodynamic limit. A different type of nonlocality 
is explored in Sec. |Vl 

Both local-density approximations, (|20p and (PT|) . as 
well as the gradient expansion (fT5|) - P^ . are applicable 
to very general classes of model Hamiltonians. In the 
remainder of this paper, we employ concepts arising in 
the context of the renormalization group, in order to ex- 
emplify and test these approximations. The concept of 
universality classes of models near criticality allows us 
to consider one representative of critical quantum spin 
chains, and deduce from it results expected to be valid for 
the entire class. Specifically, we use the LDAs (PH)) and 
pT|) to obtain the energies of the XXZ chain and test the 
possibility of obtaining the conformal anomaly and criti- 
cal exponents from LDA. The conformal invariance of the 
model allows us to obtain predictions for these quantities 
from energies of finite chains. (Later on, in Sec. |Vl con- 
formal invariance is used not only to extract information 
from functionals, but also to construct nonlocal approx- 
imations.) For Kohn-Sham Hamiltonians that belong to 



the same universality class as the interacting Hamiltoni- 
ans, universality guarantees that LDA critical exponents 
are given correctly. Nonuniversal quantities, such as the 
sound velocity, are predicted wrongly by some forms of 
LDA, but a suitable rescaling of the Hamiltonian can be 
used to cure this defect. 



III. THE XXZ QUANTUM CHAIN 

Critical quantum chains belonging to the same univer- 
sality class have the same critical exponents and, in the 
case of conformally invariant systems, the same confor- 
mal anomaly c. We expect that density functionals for 
such systems will share the same general features. In 
this section we consider the universality class comprising 
the c = 1 conformally invariant quantum chains. These 
quantum chains exhibit a critical line with continuously 
varying critical exponents. The underlying field theory 
describing the long-distance physics of such models is the 
Gaussian model (or Coulomb gas) This field theory is 
described by operators composed of a spin- wave ex- 
citation index n and a "vortex" excitation of vorticity m. 
The anomalous dimensions of these operators, related to 
the critical exponents of the quantum chain, are given by 

Xn,m^n'^x + —, n, m = 0, ±1,±2, (22) 

The parameter x varies continuously along the criti- 
cal line and its value depends on the interactions en- 
tering the particular modeli^l Examples of models ex- 
hibiting such critical behavior are the anisotropic spin- 
S Heisenberg model,— the Ashkin- Teller modeLi| the 
anisotropic triplet spin model introduced in Ref. [2^, the 
ferromagnetic spin-S' Babujian-Taktajan models,— the 
anisotropic biquadratic spin 1 chain,^^ etc. 

The most studied model in this family of c = 1 models 
is the anisotropic spin-^ Heisenberg chain, also known as 
the XXZ chain, whose Hamiltonian is given by 

^xxz(^) = -iX:K<^J+i+<^J<+i+Aa;a|+i). (23) 

Here a = {a^ ,(7^ ,a^) are Pauli matrices, L is the lattice 
size and A is the anisotropy, which in the DFT framework 
plays the role of a spin-spin interaction. This Hamilto- 
nian is exactly integrable through the Bethe ansatz.^- 
Its exact solution shows that, in the bulk limit L — )■ oo, 
the model is critical (gapless) for anisotropics — 1 < A = 
— cos 7 < 1 . As a consequence of the conformal invari- 
ance of the model, the critical exponents are exactly given 
by Eq with the model dependent parameter 

X = XA given by 

XA = - — A = -cos7, 0<7<7r. (24) 
27r 

The XXZ quantum chain is then a natural candidate to 
develop approximate functionals for DFT applications to 
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critical models in the c = 1 universality class. The ex- 
actly known finite-size corrections of the modeliS will al- 
low us to test the general ideas presented in Sec. HIl 

As in Sec. IIIBl we consider a non integrable gen- 
eralization of (|23p where besides the intersite spin-spin 
interactions we also include the on-site interactions of 
the spins with an inhomogeneous site-dependent exter- 



nal magnetic field /i^ 



1, . . . , L), namely, 



i/^^''(A, {h,}) = (^) + y 



(25) 



where we added a convenient constant. In this paper we 
study the Hamiltonian (pS)) with several types of bound- 
ary conditions. For periodic {p — 1) and open boundary 
conditions (p = 0) we have 



CTL+i =pai, p = 0, 1, 
whereas for twisted boundary conditions 



1 = e 



'L + l 



(26) 



(27) 



where < </> < 27r and ct^ — \{(t^ ± icr^) arc the usual 
raising and lowering SU(2) spin operators. 

The inhomogeneous Hamiltonian (|25p is suitable for 
the DFT approach presented in Sec. |TT1 It has a C/(l) 
symmetry, and it is easy to obtain a simple Hamiltonian 
H", as in ([7]), with the properties (i)-(iii) discussed in 
Sec. ini The Hamiltonian (|25p commutes with the global 
charge 



Hs are given by arbitrary combinations of free fermion 
energies, given in terms of the eigenvalues of a L x L 
matrix, whose elements depend on the boundary con- 
dition in ([50]) and the values of the magnetic fields /if 
{i = 1, . . . , L). This implies that with reasonable com- 
puting efforts we can diagonalize Hg for lattice sizes up 
to L = Lniax ~ 5000. This is certainly much more we 
could reach for the 2^ x 2^ interacting Hamiltonian (|25p . 

The inhomogeneous magnetic field hf [i — 1, . . . , L) 
are fixed, as in (jl2p . by imposing that the ground state 
eigenfunction of ([50)1 and (pS)) . on a given sector with 
fixed number N of up spins, share the same density dis- 
tribution rii (i = 1, . . . , L) of up spins: 



dW\ 



dn 



i = l. 



(31) 



As in ((33)) . W[n] is given by the difference of the function- 
al of the interacting (O and non-interacting ^ models. 
Since in the present case Hg is obtained from H by set- 
ting A = 0, the corresponding functionals are related and 
we can write 



W[n] = FA[n] - Fa^qM, 

where Fa is the functional obtained by using in 
Hamiltonian iJ^^^(A) given in i.e., 

FaM = min|*>„^ „^ < ^\H^^^{A)\^ > . 



(32) 



the 



(33) 



The density profile nP = (n^, 



of the ground state 



of H-^^'^{A, {hi}) restricted to the sector with N up 
spins, is obtained by solving self-consistently ((30| with 
([?T|) . From (fT51) the ground state energy of the inhomo- 
geneous interacting model (^5)1 is given by 



. . af + l . ^ 
rii. Hi = — - — , I = 1 



(28) 



1=1 



which gives the total number of up spins in the cr^-basis. 
By comparing ([1]), ([2]) and ([25]) we identify 



ycxt ^ ^ 



-h,, i^l,....,L. (29) 



The simple Hamiltonian that plays the role of the Kohn- 
Sham Hamiltonian is obtained by setting A = in (pS)) . 
i.e., 



En = E[n°] = Esiiht}) 



L 

E 

i=l 



^-n° + W[n°]. (34) 



Here Es{{hf}) is the ground state energy of the free 
fermion XY Hamiltonian H^^^{A = 0, {/if}), with N = 

Z]f=i '^i spins. 

All difficulties in obtaining the energies En in 
now reduce to the derivation of the functional Wa [n] in 
Following the approximations considered in Sec. [Tll 
we next consider the XXZ chain with distinct boundary 
conditions and inhomogeneities. 



H, {{ht}) = i/^^2(A = 0, {/if}) = i/^^^(0) - E 

(30) 

This simple Hamiltonian is the well known XY model in 
the presence of site-dependent magnetic fields /if (i — 

I, . . . , L). Through a Jordan- Wigner transformation^^ 
Hs is tranformed into a Hamiltonian describing L non- 
interacting spinless fermions. The 2^ eigenenergies of 



IV. LOCAL APPROXIMATIONS FOR THE XXZ 
CHAIN 

In the absence of external fields /i^, the only source 
of inhomogeneity (i.e., breaking of translational invari- 
ance) are the boundaries. In this section we consider the 
Hamiltonian (|23p on finite chains, for periodic, twisted 
and open boundary conditions. 
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A. Periodic boundary conditions 

In this case, beyond the conservation of the total num- 
ber of up spins, the Hamiltonian is also trans- 
lationally invariant. As a consequence, we can sepa- 
rate the Hilbert space associated with the interacting 
model H^^'^{A), as well as that of the auxiliary Kohn- 
Sham Hamiltonian H-^^^{A = 0), into block-disjoint 
sectors labelled by the density of up spins n — jj, 



(n — 0, -J-, . . . , 1) and momentum P — 77, 

0. 1....,L — 1). In each of these sectors we can apply 
the DFT-LDA procedure of section [ll] to calculate the 
lowest energy Eo{L,n,p). 

Let {(j) > he any vector on the arbitrary sector {n,p), 

1. e., J2i=i f>-ilL\4> >— M't' > T\(f> >— e'^P\(f> >, where 
T is the translation operator. It is simple to show that 



<(f)\hi\(f)> J2i 



< 



> 



L 



N 
1 



(35) 

This implies that any eigenvector of H^^^(A) has a uni- 
form distribution of densities = n (i = 1, . . . , L). 

The functional i^A = FA{n,p), for general values 
of A, is in the present case a simple function of the two 
variables n — N/L and p. Consequently the local mag- 
netic fields {hj} in the auxiliary Hamiltonian ((30)) are 
site independent: 



dn 



dn 



{FA{n,p)-Fo{n,p)) = h' 



(36) 

Hs{{hf}), in this case, is just the XY model in the pres- 
ence of a uniform magnetic field hg . The lowest energy of 
this auxiliary Hamiltonian, in the sector {n,p), is given 
by 



Esih") = Le"A^S{n,p) - Lh^n, 



(37) 



where e.^^^ iji^p) = Cq '''°'^(n,p) is the lowest eigenenergy 
per site in the sector {n,p) of the XY model in the absence 
of external magnetic fields. 

In order to obtain the lowest eigenenergies of the in- 
teracting model on the sector {n,p), from (p4)) . it is nec- 
essary to know the functional WA{n,p). Following the 
discussions of Sec. IIIIl we use the local approximations 
(PO)) and (HD) for this functional. In the present case they 
are given by 



2^1 



L,pcr / \ r L,pcv 

A [n,P)^Lw^'' 



{n,p), (38) 



W^''^'^{n,p) = ^u;^(n,p) = Lw^{n,p). (39) 



1=1 



Here w^^"^" {n,p) 



— e 



L.pcr 



{n,p) 



L,pcr 



{n,p) and 



^ {n,p) = e'^{n,p) — eg°(n,p) are the interaction ener 
gies per site, and e^^'^'^{n,p) and eA{n,p) represent the 



lowest total energy per site in the sector {n,p) of the pe- 
riodic Hamiltonian H^-^'^{A,p) with L sites and with an 
infinite number of sites, respectively. The Hamiltonian 
i7^^^(A,0), with periodic boundaries, is exactly inte- 
grable through the Bcthc ansatz.^^ The eigenenergies for 
lattice size L, are obtained by solving a set of = nL 
coupled non-linear equations. This can be done,^^ at 
least for the lower eigenenergies of sectors {n,p), either 
for lattice sizes up to L « 1000 or in the bulk limit 
{L 00). 

Suppose, for the moment, that we have evaluated ex- 
actly e^^°''{n,p) and the LDA functional ([38]) and p9|l . 
The lowest eigenenergy, in the sector (n,p), is obtained 
by inserting (j37p in (|34p and using one of the approxima- 
tions ([55]) or ([5^ . From the finite-size LDA fimctional 
(l38l) we obtain 



E^^'^ {L,n,p) = Lcq'^'^'' {n,p) — Lh'^n 



^ dWA{n,p) 
y — 1- WAKn.p). 



dn 



(40) 



The relation ([5^ gives, in this case, the exact result 



Eo^^'^{L,n,p) 



Le 



L.pcr 



{n,p). 



(41) 



This result is expected, since the density distribution of 
up spins for any of the eigenlevels is spatially homoge- 
neous, and we only need to keep the first term in the 
expansion (|15p - p^ . The conformal anomaly and criti- 
cal exponents that are calculated from the leading finite- 
size corrections of the eigenenergies, are then obtained 
exactly. 

This is not the case for the LDA functional ((39)) , which 
derives from the infinite chain. If we use ([55)1 and (|57)) 
in (|M)) . we obtain 

E^^'°° {L,n,p) — LeQ'^'^\n,p) — Lh^n 



|][— (e^(n,p)-eg°(n,p))]n 
+L{e^{n,p)-e^{n,p)). 



(42) 



Equation (|36)) gives 



L[e^{n,p) 



E 

L,pcr 



LDA. 



\L,n,p) 



(n,p) -eg°(n,p)]. 



(43) 



This result is certainly not exact for all orders of 1/i, 
unlike that obtained from the finite-size LDA functional 
psp . However, in typical applications of LDA in ah ini- 
tio electronic-structure theory, as well as to the Hubbard 
and the Heisenberg model, one employs an LDA func- 
tional deriving from the corresponding infinite system, 
for any calculation on finite-size systems. In the present 
work, this approach is represented by the infinite-size 
LDA functional (|55)) . It is therefore interesting to eval- 
uate the leading finite-size corrections of (|^5)) . obtained 
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from (|39p. to see if we can, from this more standard LDA, 
still obtain exact results for the conformal anomaly and 
critical exponents. 

The ground state of the XXZ chain has zero momen- 
tum. It belongs to the half- filled sector with n — — ^, 
for even values of L, and to the sector n = ^ ± in the 
case of odd values of L. Let us restrict ourselves to the 
case where L is even. The asymptotic behavior of the 
ground state energy is exactly knowni^i^S 

i?r(^;^,0) = Lei'P-(i 0) 



where ca = 1 is the conformal anomaly and 



TT sm7 

7 



-1< A 



' cos 7 < 1, 



(44) 



(45) 



is the A-dependent sound velocity. Using (j44|) with A = 
in we obtain the LDA prediction for the asymptotic 
behavior of the ground state energy 



E 



LDA.oo 



(i;i 0) = Le^(i 0)-^CAV6i + o(|), (46) 



where ?;o = wa=o — 2. 

Upon comparing (|46p with (|44ll . we see that the LDA 
([55)1 yields an incorrect amplitude for the 0{1/L) term. 
The sound velocity predicted is not given by (|45|) but by 
its value at A = (7 = ^). This means that the LDA 
([55)1 does not allow an exact prediction of the conformal 
anomaly. 

The critical exponents, as a consequence of the con- 
formal invariance of the infinite system, are obtained 
from the finite-size corrections of the energy gaps of the 
quantum chain. The lowest eigenenergies on the sectors 
with = "I + i/ (i^ = 0, ±1, ±2, . . .) up spins (density 
n = = ^ + ■^) are zero- momentum states. The asymp- 
totic behavior of these energies is exactly knowni^i^ 



i?A(i;;^ + T>0)-i4'P"(- + ^,0) = 



1 V 



— VAXA'^^ + 0{-), (47) 



where xa is given by ([24| . The energy gaps related to 
the eigenenergies (|T7)) then give us the parameter xa that 
characterizes the particular c ~ 1 conformal theory: 

G'i = EAiL; i + 0)-EAiL; i, 0) = "^VAXAi^^+oi^), 

(48) 

where 1/ = ±1 ± 2, . . .. 

The gaps (gH]) predicted by the LDA functional ^ 
are obtained from 

- 



L.pcr 



2' 

+i(e^(^ + ^,0)-eJ(i,0)) 
-i(eS°(^ + ^,0)-eo-(i,0)). 



(49) 



The leading behavior, as L ^ cxo, of the energies in the 
first parenthesis of the last expression is given by ([T7|) 
with A = 0. In order to find the leading finite-size cor- 
rections of the gaps (j49|) . we need to calculate those cor- 
rections for e^(i -1-^,0). This is the lowest eigenenergy 
of the infinite system with finite density n = i + of 
up spins and zero momentum. It is important to stress 
that i^/L is small but finite. From (|47| we have the exact 
behavior 

eA(UT'0)= lim £;a(L', 1 + ^, 0)/i' 



27r ,L'v 



1 



(o'O) + —^i — yXAVA + o{-^ 



,1 



e?(:.,0) 



By substituting ([501) and (gZl) in 



^TTVAXA 2 , /In /rn\ 

^ +o(-^). (50) 
we find 



G 



LDA,i 



2 . 1 N 

—VAXAf 



±1,±2, 



(51) 



which reproduces the exact result This means that 

although the infinite-size LDA ([39)1 does not give an ex- 
act value for the conformal anomaly c, it produces exact 
results for the critical exponents xav^ ■ In order to ob- 
tain the exact results (|?T|) . it is crucial to use the exact 
result for the energy per site of the homogeneous model, 
e^(n,p), which can be obtained from the Bethe ansatz 
solution of the model. For completeness we give, in 
the appendix, the relevant integral equations that pro- 
duce exact values for e^(n,p). In Sec. |V]we introduce 
an analytical approximation for the functional (j39p that 
contains all relevant ingredients to reproduce the exact 
results of the critical exponents given in ([5T|l . 

Let us now consider the LDA predictions for the mass 
gap amplitudes related to the sectors with momentum p. 
The energy-momentum dispersion relation gives 

Gl = L(ei-P"(i; ip) - el ^"(L; 1,0)) = vaV. (52) 

where va is the sound velocity (|45p and V — ^ 
0,1,. 



From this last expression we also obtain 



1 



e^(^,P)-e^(^,0) 
1 



lim -1[L'(4' P-(L'; ip) - ef' ^^'■(L'; i, 0))] 

L'-^oQ L/ Z Z 

= lim VAP/L' = 0. (53) 

L' — >oo 

From (gSl), ([52]) and ([53]) we obtain the LDA prediction 
for the momentum mass gap 



G 



LDA.p 



(54) 
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Again, as in (|46|) . we obtain the same incorrect predic- 
tion for the sound velocity of the model, i.e., the value of 
the non interacting XY model. However, the sound ve- 
locity is a model-dependent quantity. Its value changes 
if the Hamiltonian is multiplied by an arbitrary posi- 
tive constant, without changing the ratios of mass gaps 
that are related to the critical exponents. On the other 
hand, the underlying conformal field theory governing 
the critical line of the model should have a fixed sound 
velocity. This expectation is bourn out by the XXZ quan- 
tum chain if we consider in (|23p the rescaled Hamiltonian 

H {A)/vA, where va is given by ([45|1 . In this case 
the sound velocity becomes unity for any value of A and 
we obtain, even from the infinite-size LDA (j39p . exact re- 
sults for the sound velocity {va = 1), conformal anomaly 
(c = 1) and critical exponents (xa^^^). 

Before finishing our discussion of the periodic case, it 
is also interesting to consider the mass-gap amplitudes 
related to the special sector with momentum p = tt. The 
numerical solution of the Bcthc ansatz equations for fi- 
nite chainSfi^ and the conformal invariance of the infinite 
system predict, from the finite-size corrections related to 
this gap, the critical exponent xq^i = 1/4:Xa, where xa 
is given by ([M)) . As before, the infinite-size LDA ([5^ 
reproduces this exact result. 



B. Twisted boundary conditions 

In this case, the boundary condition (|27|) is specified 
by the angle < (j) < 2tt, the periodic case correspond- 
ing to (j) — 0. The total density n of up spins remains 
a good quantum number. Moreover, the quantum chain 
also displays a generalized translation invariance for ar- 
bitrary values of (p (see Ref. 'sl] for the proper definition 
of translations on the lattice) . As a consequence of this 
invariance, as in Eq. (|45)) . the density obtained from any 
given eigenstate of the Hamiltonian is homogeneous, i.e., 
rii = N/L = n {i = 1, . . . , L). For the sake of simplicity, 
in the following we restrict ourselves to states with zero 
generalized momentum. The LDA functionals, obtained 
from and (PU)) . are generahzations of ([55]) and 
and now take the form 

ly^D^'^n) = Lw^^'^in), (55) 

and 

W^A°^'°°(^) = i^«A(«)- (56) 

Here w^^{n) = e^'^{n) — e^'''^{n) and w^{n) = e^(n) — 
ej^(n) are the interaction energies, per site, in the eigen- 
sector with density n of up spins and lowest energy 
e^''^(n) or e'^{n) of the XXZ chain with twisted bound- 
ary conditions specified by the angle 4>. 

Let us restrict ourselves to the case where L is even. 
The ground state energy belongs to the eigensector la- 
belled by the density n = 1/2 of up spins. The Bethe 



ansatz gives the asymptotic behavior of the ground state 
energy for L ^ oo, 

Ei{L, 1) = Le^{\) - ^~CA + o(i). (57) 

Here 6^(5) is the ground state energy, per site, in the 
bulk limit, and 

6a = 1 - 122;„ i = 1 - 12(^)V4a;A, (58) 

• T TT 

is the conformal anomaly of an effective model described 
by the XXZ quantum chain with twisted boundary con- 
dition specified by (f>M- The parameters Xn^m, va and xa 
in Eqs. ^ are given by Eqs. ([221) , dill) and dH]), re- 
spectively. The exact resultsi^ for the leading behavior of 
the mass gaps of the eigensectors with density n = ^ + -j- 
(i^ — ±1, ±2, . . .) of up spins yield 

= EiiL, \ + ^)- EiiL, i) ^ ^v^xAi^^+oi^l 

(59) 

independently of the boundary angle (j). 

As in the periodic case, we now test the possibility to 
recover the asymptotic behavior specified by Eqs. (|57p 
and from the LDA functionals ([HS]) and ([51]). Simi- 
larly to the case = (periodic boundary conditions), if 
we use the finite-size LDA ((55|) we obtain the exact result 

i?5'°^^^'°°(i,«) = iei'^n), (60) 

and thus the correct value for the effective conformal 
anomaly ca and the critical exponent xav^ ■ This hap- 
pens because in this case ([55)) is exact. For any eigen- 
state of the Hamiltonian the density is homogeneous 
(n\ ~ ■ ■ ■ = til) and the only nonvanishing term in the 
expansion (|17p is the first one. 

Let us now consider the infinite-size LDA functional 
([55)) . The expressions for the LDA eigenenergies, which 
for the periodic case (0 = 0) were given by ([43)) . are now 
replaced by 

ijLDA,oo(^) = L[et{n) -f e^-^n) - e-(n)]. (61) 

By evaluating Eq. ([57)) for A = we obtain from Eq. ([IT)) 
the asymptotic behavior of the ground state energy, as 

L ^ 00, 

^LDA.oo( 1 ^ = ie^^^^, 0) - ir-CAVnl&L + o(l/L). (62) 

Comparing ([62)) with ([57)) we obtain, as in the (j) — Q 
case, the sound velocity iiq = 2 of the XY model instead 
of the correct VA^ given in ()45p . 

As in the periodic case, in order to evaluate the critical 
exponents, we need the asymptotic large-L behavior of 
e^(i + f), with f small but finite {v = ±1,±2,...). 
This is obtained from ()59p in a similar way as we got 
([50| from ([47]): 

+ j) = + 2t,vaxai^^IL^ + o(1/l2). (63) 
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By inserting the LDA energies ((611) in ((59)) . and using 
we obtain 

^LDA,..oo ^ ^„^^^^2 ^ (g4) 
ij 

Thus, as in the periodic case, the infinite-size LDA ((56|) 
gives exact results for the critical exponents, although it 
is not the exact functional for the finite-size chain. 

The wrong estimate obtained for the sound velocity is a 
consequence of the distinct sound velocities ua and uq — 
2 of the interacting and non-interacting Hamiltonians. If 
we consider the rescaled Hamiltonian H-^'^'^{A)/va, sim- 
ilar to what we did in the periodic case, the sound velocity 
is now unity for both Hamiltonians, and the prediction 
from the LDA is now exact both for the conformal 
anomaly and critical exponents. 

C. Open boundary conditions 

For open boundary conditions, the quantum chain (1231) 
still conserves the total number TV of up spins, but now 
is not translationally invariant. In general, the eigen- 
functions belonging to the eigensectors with total density 
n = N/L of up spins produce inhomogencous distribu- 
tions of local densities ni, . . . , n^. Let us restrict our- 
selves, in the following, to cases where the lattice size L 
is an even number. 

The Hamiltonian (|23p with open boundaries commutes 
with the spin reversal and total z-magnetization opera- 
tors, 

L L 

R = l[<7f, S^ = Y.^t- (65) 

i=l 1=1 

When restricted to the sector where = 0, these opera- 
tors also commute with each other, [R, S] = 0. This im- 
plies that in the eigensector with density n = N/L =1/2 
of up spins, the parity under spin reversal is also a good 
quantum number. It is then simple to show that an ar- 
bitrary eigenfunction >, on this sector, produces a 
homogeneous local density, namely, 

=< + > / < ^-1^- >= i, z = 1, . . . , L. 

(66) 

Here we will only consider this case, whose analysis is 
similar to that presented for periodic and twisted bound- 
ary conditions. The general case, in which the density 
is spatially inhomogencous, is treated by means of an 
improved functional in Sec. IV Bl 

The ground state of ((23|) belongs to the special sector 
where n — N/L = 1/2. The XXZ quantum chain with 
open boundaries is also exactly integrable through the 
Bethe ansatz.'^? The exact asymptotic behavior of the 
ground state energy is known'^^ to be 



= Lei'°P- = Le^(i) + /A-^^+o(^), (67) 

where ca = 1, is given by (|i5)) . and is the surface 
energy due to the open ends of the lattice. According to 
conformal invariance^ the leading behavior of the mass 
gap associated to the eigensector with density n = ^ + ^ 
(y = ±1, ±2, . . .) of up spins is given by 

= -^wasja" + o(-), (68) 

where a;^" is a surface critical exponent. The finite-size 
analysis of the model^S. gives the exact result 

x'^^ = 2xisj/^ , (69) 

where xa is given by ((24|) . Comparing (|68p with (jSQ)) 
and using (|69p we verify that the lower mass gaps {y — 
±1, ±2, . . .) does not depend on the boundary condition, 
up to order (1/i). 

We next verify whether the exact results ([ST]) and ([SS]) - 
([69)) can be obtained from the LDA functionals intro- 
duced in Sec. Illli whose infinite-size and finite-size ver- 
sions now become 

L 

W LDA, L, open/ \ \ ^ L,opcn/ \ incW 

^ ''^ (m, . . . ,nL) = 2^ w^' ^ (rij), (70) 

1=1 

and 

L 

I^^°^'°°(ni,...,ni)EE^w;^(nO. (71) 

i=l 

When the latter approximation is applied to an open sys- 
tem, it misses corrections of order L", which arise from 
the fact that in a finite-size open system one nearest- 
neighbor interaction is missing, as compared to the finite 
periodic case. In the thermodynamic limit the resulting 
error is negligible, but for finite systems it can be approx- 
imately corrected by replacing the preceding equation by 

r - 1 ^ 

<°^'°°(ni,...,ni)«-^^w^(n,), (72) 

4=1 

which we take as our definition of W^^'"^ in this case. 

The ground state belongs to the special eigensector 
n — 1/2, where we have a homogeneous distribution of 
densities of up spins. This means that, as in the periodic 
case, use the finite-size LDA ([7Dt yields the exact result 

^LDA,L,opcn(^^ 1) = Lei'°P™(i), (73) 

implying exact predictions for the conformal anomaly, 
sound velocity and surface energy. It is interesting to 
mention that in order to obtain the exact result (I73p it 
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was crucial to use e^'°''™(n) in the LDA ([70)) . If we used 
(L — l)e^^°'^ / L instead, we would get an inexact result, 
and wrong predictions for the sound velocity and surface 
energy. 

Let us now consider the results for the ground state 
energy following from the LDA ([7T|) . In this case we 
can derive the results from those obtained from the LDA 
(I39p . appropriate to the periodic case. The external mag- 
netic field (f3T|) used for the Kohn-Sham auxiliary Hamil- 
tonian (|30p is the XY model with open boundaries in the 
presence of an uniform magnetic field. The ground state 
energy of the Kohn-Sham Hamiltonian is now given by 

EM = Le^AT\l)-Lh^l- (74) 
By plugging the LDA ^ in (j34p we obtain 

+ ^(^A (^) - ecTC^))], (75) 

where we have used (j74p and the magnetic field given in 
([55)1 . The asymptotic behavior (|^ and (pT]) gives us, as 

L CXD, 

+ + + (76) 

Upon comparing ([7S)l with the exact result (|F7|) . we see 
that the LDA ([7T|) gives wrong results for the surface 
energy, sound velocity and conformal anomaly. 

The sound velocity va and the surface energy are 
non universal quantities. Consequently is not a surprise, 
due to our experience with the periodic case, that the 
LDA ((7T|) does not give exact predictions in the present 
case. The underlying conformal field theory governing 
the fluctuations of the quantum chain is defined on a 
half plane and has a fixed sound velocity and surface 
energy.i^^^ We can also obtain a quantum chain with 
this property by considering the rescaled Hamiltonian 

J?^X2(A) = (i/^X2(A) - Le^ - fD/vA. (77) 

where H^^'^{/S.) is the Hamiltonian with open 

boundary conditions. For this Hamiltonian, the ground 
state energy e°° = 0, the sound velocity va = ^ and the 
surface energy = are constants. The infinite-size 
LDA ([7T|) . applied to (|77p. recovers the correct results. 

In order to evaluate the energy gaps, in the case of open 
boundaries, we need to consider the lowest eigenenergy 
in the sector with average density ^ — ¥^ In this 
case the density distribution is not homogeneous anymore 
and we must diagonalize the Kohn-Sham Hamiltonian 
numerically, and solve self-consistently for the densities. 
These calculations are presented in Sec. IV Bl Before 
performing these calculations, we introduce, in Sec. IV Al 
approximate analytical expressions of the functionals. 



V. NONLOCAL APPROXIMATIONS FOR THE 
XXZ CHAIN 

In the preceding section we applied the finite-size and 
the infinite-size LDA functionals to systems in which the 
only source of inhomogeneity are boundaries. The finite- 
size LDA, due to its nonlocal dependence on the system 
size, turns out to be superior, but even the more conven- 
tional infinite-size LDA produces correct critical expo- 
nents, and, for rescaled Hamiltonians, also correct con- 
formal anomalies. However, neither LDA is guaranteed 
to be reliable for systems in which the density distribu- 
tions is inhomogeneous also in the bulk, as is the case in 
the presence of site-dependent external fields h^^^. We 
now develop and test nonlocal approximations for such 
inhomogeneous systems. 



A. Construction of nonlocal functionals 

In order to calculate the ground state energies pro- 
duced by inhomogeneous density distributions it is nec- 
essary to solve selfconsistently the Kohn-Sham Hamilto- 
nian dSOD with /if given by ([21])- ([22]). The quality of the 
results will depend on the approximation used for the 
functional W, i.e., on the number of terms kept in the 
expansion (|15p. 

In section ITVl we considered only the first term in (|15p . 
producing the functionals ([381)- (011), (EUl-lIMl) and ((TOl) - 
(I7ip . for periodic, twisted and open boundary conditions, 
respectively. These functionals are exact for completely 
homogeneous systems, and produce reasonable approx- 
imations for weakly inhomogeneous situations. To im- 
prove the results also in more strongly inhomogeneous 
cases, we now additionally consider the second term in 
(fTS)) . As we shall see, this extra term is necessary for the 
correct evaluation of the mass gaps of the quantum chain 
with open boundaries. 

On periodic lattices, the functionals exhibit general- 
ized translation invariancc in the sense that for any inte- 
ger TO, Fa{ni, . . . , ul) = Fa{ni + m, . . . ,nL + to), where 
a — I, s, 7ij+L = nj (j = 1, . . . , L). In this case the func- 
tions ai[n] ^ ■ ■ ■ = aL[n] ~ a[n] defined in are site 
independent and the second term in (|17p reduces to 

A TV 

^aj[nj]{nj - Ui) = ia[n,;](— - n^). (78) 

From P^ -(|19 p we then have 

iy^°^^°°(ni,...,ni)= (79) 

^ ^ N 

^zi;^(n)|„=„, +^ia[n,](--nO, (80) 

i=l i=l 

where X^f^i ni = w'^{n) = e^(n) - eg°, e^(n) is the 
ground state energy per site of the infinite homogeneous 
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system with total density n of up spins, and 

1 H/a(- ■ • , "-i, ni + 5, rii, . . .) - M^A(ni, . . . , rii) 



afnil = lim 

s^o L 



In order to estimate a[n] we approximate 



(81) 



WA{---,ni,ni+6,ni,...) w (L-l)w^(ni)+w^(n,+7<5/2), 

(82) 

where 7 « 1. The last term gives the contribution due 
to the inhomogeneity ni + 5. Since the Hamiltonian is 
composed of two-body interactions we expect that the 
contribution due to the inhomogeneity is given by the 
energy evaluated around the average density (n^ + (n,; + 
(5))/2, and hence^ "f k, 1. By plugging 
obtain 



n we 



a[n,\ = -- 



dS 



2L drii 



(83) 



and hence, 



1=1 



TjrLADA.oo / 



9n 



(") 



(84) 



In the case of open boundaries, differently from the 
periodic case, the number of links on the lattice is (L — 
1). Moreover, the second term on the expansion in (|17p 
is now site dependent. An approximate functional for 
this case is obtained by taking the contributions of the 
boundary sites 1 and L to be half that of the other sites 
i — 2,...,L — 1. Similar arguments as those used to 
obtain ([7T|) for the first term and ([5^ for the second 
term in (jl7p . then give the functional 



(ni,. 



,"-l) = 



L-1 
L 



dw°. 



dn 



■(n) 



(85) 



Approximations (I84p and ([85)) . for periodic and open 
chains respectively, can be used, up to this point, for any 
quantum chain. Note that although we started out from 
a functional depending on the lattice gradient n,j —n.i, the 
final expressions depend only on the differences N/L — rii, 
i.e., on the deviation of the local density from the average 
density of the system. This average makes these function- 
als highly nonlocal, in a similar way as described below 
Eq. (|2ip for the finite-size LDA. The present function- 
al , however, depend on L in a more complex way than 
the finite-size LDA, and additionally depend on N. For 
these reasons we refer to ([M)) and ([55]) as local and av- 
erage density approximation (LAD A), highlighting thus 
its simultaneous local and nonlocal dependence on rii. 
(In fact, there is a remote conceptual similarity to the 
average-density approximation of ab initio DFT,^ which 



also displays a nonlocal dependence on the density via 
an integral of n(r) over a certain region in space.) 

To make these functionals useful in practice, it is nec- 
essary to evaluate the ground state energy per site (n) 
for fixed density n. Actually this gives us a general pro- 
cedure to produce a good approximate functional for an 
arbitrary non integrable quantum chain. Numerical di- 
agonalization of small chains, through the Lanczos or the 
density matrix renormalization group (DMRG) methods, 
can give us an estimate for the ground state energy as a 
function of the density. This estimate is then used to pro- 
duce the functional ([M)) and 6.6. The advantage of study- 
ing a chain whose homogeneous limit is exactly integrable 
is precisely the exact evaluation of these quantities. For 
the XXZ chain, they are obtained by solving the coupled 
integral equations (|A.ip - (|A.3p of the appendix. We could 
solve these equations numerically for a set of discrete den- 
sities and use the obtained results to numerically define 
e^(n) for arbitrary densities < n < 1. Alternatively, 
we can use all known information from the exact solu- 
tion and the conformal invariance of the XXZ model, to 
produce a good analytical parametrization, e^'^{n), for 
e^(n). Below we follow the second, analytical, route, 
but test the resulting parametrization by comparing it 
to numerical data. 

The Bethe ansatz solution of the XXZ chain at A = 
and A — > — 00 indicates that at those couplings the model 
describes spinless non interacting fermions (up spins). 
The model with A = is the standard XY model where 
the up spins have a hard-core interaction of unit size, 
in lattice space units, that forbids double occupancy of 
up spins on one lattice site. The model with A —>■ —00, 
on the other hand, is related to an effective model that 
forbids the occupation of pairs of up spins at distances 
smaller or equal to unity, in lattice space units 1^ Effec- 
tively, apart from harmless constants, these models de- 
scribe spinless particles with hard-core sizes s (s = 1 for 
A = and s = 2 for A —00). In these cases an ana- 
lytical solution, for arbitrary densities < n < 1, can be 
derived^^ 



-(1 - (s - l)n) sin( 



where s = 1 for A = and s — 2 for A — > — 00. On the 
other hand, for —1 < A < 1 the model is gapless and 
the ground state belongs to the sector with = L/2 up 
spins, or average density n = N/L = 1/2 of up spins. 
Spin-reversal symmetry gives 



':{l-n), 0<n<l. 



(87) 



Moreover, for arbitrary values of — 1 < A < 1 conformal 
invariance implies for the eigensectors with N — ^ + v 
(ji^l << L) of up spins the leading finite-size behavior 



L 
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or equivalently 



in) 



l)2+o(-l), (89) 



for n w 1/2. The parameter a; a is given by Eq. (|24p and 
gives the critical exponents (see Eq. (|22p ). and va is the 
sound velocity given by Eq. (|45p . This means that 



de 



L,per. 



n) 



0, '^'^ 



(n) 



d^r 



(90) 



CoUecting the information in Eqs. ([55]) . ([57]) and 
we obtain an approximate analytical parametrization of 
the ground state energy at arbitrary densities for — 1 < 
A < 1 



=P^'^(n) = --(l-(s-l)n)sin( 



1 - (s - l)n 



and 



A(l-An) + ci + ^C2in-^f, in<^), (91) 



^(n)=er(l-^), in>h. (92) 



TABLE L Parameters defining the parametrization e^^"^ given 
in (|9ip - (|92p for several values of the anisotropy A = — cos 7. 



7 

Sir 
6 

3ir 
4 

2tv 
3 



A 

0.866025 
0.707107 
0.5 

-0.5 
-0.707107 
-0.866025 
-1 



-1.691654 
-0.385388 
0.295398 
1 

1.349051 
1.454903 
1.527244 
1.583761 



C2 

0.326211 
0.444007 
0.389751 


0.323765 
1.133043 
2.395693 
5.625124 



Cl 

-0.013096 
-0.013299 
-0.009221 


-0.003210 
-0.010000 
-0.017932 
-0.026728 



Fig. [2] ("-) and e^^'^(n) at anisotropics A 



and 

A = We see from these figures that for 0.4 < 71 < 0.6 
the largest deviation is around 0.2%. This means that for 
inhomogeneous density distributions {ni,i — 1,...,L} 
with 0.4 < fii < 0.6 the exact result e^in) can be re- 
placed for our parametrization e^^"^ of the functionals 
(|84l) and (|85|l . with excellent accuracy. For large inho- 
mogeneities one must numerically solve the Bethe ansatz 
equations, requiring additional computational effort. 



The parameter s is obtained by imposing the first condi- 
tion in (1201) J i-e., by solving the equation 



1 . , TT 

-sin(- 



:) 



■cos(: 



A = 0. 



(93) 



The parameter C2 is obtained from the second condition 
in (190]) . namely 



C2 = AnVAXA 



167r TT 

■ sm( 



(3-.)3 



^3 - s' 



(94) 



Finally, the parameter ci (a harmless constant for the 
evaluation of gaps) can be obtained by imposing that at 
n = i, e^^'(i) coincides with the exact result e^(i). 
This is obtained by solving the Bethe ansatz equation in 
the appendix with A —^ 00. 

The parametrization (I9ip - (I5^ is then obtained from 
the solution, for each A, of ((93)) (giving the parameter 
s) and by solving the Bethe ansatz equations (|A.ip - (jA.3p 
at the density n ~ ^. This procedure is certainly much 
simpler than the brute force numerical solution, for each 
density n, of the integral equations (|A.ip - (jA.3p . The pa- 
rameter s in (PT1) - ([5^ plays the role of a generalized in- 
teraction range, or effective size, of the hard-core inter- 
actions among the up spins. It varies from s = 1 to 
s = 2 as A goes from to —00. As A decreases from 
zero, this parameter is smaller than unity, reflecting the 
decrease of the repulsion among the up spins, since now 
the static term, controlled by A, is attractive. In table U 
we show, for some values of A, the parameters entering 
the parametrization (PT|) - ([5^ . 

In order to compare qualitatively our parametrization 
(PT|l - (|92p with the exact result we show in Fig. [T] and 




FIG. 1: (Color online) Parametrization of (n), given by 
Eqs. (|9ip - (|92p . compared to the exact ground state energy, 
e2'(w), for anisotropy A — —1/2. 

It is important to stress that the density of the interact- 
ing system is reproduced via the site-dependent magnetic 
fields {/if} of the non-interacting Kohn-Sham Hamilto- 
nian ([50)1 . These effective fields depend on the difference 
of derivatives of the functionals at the anisotropy A and 
A = 0, i.e., 



^ TjrLDA.OO/ \| 



(95) 



In order to compare the effective magnetic fields pro- 
duced by our parametrization e^^(n) with those ob- 
tained from the exact values for e^(n), we show in 
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FIG. 2: (Color online) Parametrization of e^^{n), given by 
Eqs. (|9H) - (|92|) . compared to the exact ground state energy 
eA («) for anisotropy A = 1/2. 

Fig. [3] these fields, for /i™* = 0, for two values of the 
anisotropics, in a system with constant density ni — 
■ ■ ■ = riL = n. In this case hi, ... = . The exact 
results are obtained within computer accuracy from the 
numerical solutions of the Bethe ansatz equations given 
in the appendix, and the derivatives were obtained by 
using a cubic spline fitting of the numerical data. We see 
from this figure that for 0.4 < n < 0.6 our parametriza- 
tion (PT|) - ([5^ produces a rather good approximation for 
the local potentials in the Kohn-Sham Hamiltonian. 

It is interesting to observe from Fig. [3] that the mag- 
nitude and even the sign of ft,* depend on the site den- 
sity and on the value of the anisotropy A. In Fig. 2] 
we show the effective magnetic field h'^ {n) obtained from 
the parametrization (|5D|) - (PT|) for several values of the 
anisotropy A. We see from this figure that while for 
A < the magnetic field /i* increases with the density 
n of up spins, for A > we have the opposite behavior. 
This is quite reasonable, since for A < (A > 0) the net 
effect of the anisotropy is to increase (decrease) the en- 
ergy as the density increases and the up spins are forced, 
in competition with the kinetic energy, to occupy nearest 
neighbor sites. 

Next, we apply the LADA functional dUl) and 
with ([50)1 - (PT|) . to several quantum chains whose ground 
state exhibits a spatially inhomogeneous distribution of 
up spins. Such inhomogeneities can be produced either 
by boundary conditions or by the presence of inhomoge- 
neous external fields. We are going to consider example 
of both cases in the following subsections. 

B. Open boundary conditions 

As discussed in Sec. IIVCI except for the eigensector 
with total density n = 1/2 of up spins (L even), the 



FIG. 3: (Color online) Effective field h" = h{ ^ ■ ■ ■ ^ hi 
obtained from (|95|) with = as a function of the density 
n = ni — ■ ■ ■ — riL for two values of the anisotropics. The ex- 
act and approximate curves are obtained by using in (|95|) the 
parametrization (|9ip - (|92p and the exact result, respectively. 
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FIG. 4: (Color online) Effective field h" = hi = ■ ■ ■ = hi 
obtained from (|95[) with h°^^ = as a function of the density 
n = ni — ■ ■ ■ — riL, for several values of the anisotropics. The 
curves were obtained by using the parametrization (|9ip - (|92p 
in dHSl) and ([95)). 



eigenstates of the XXZ chain with open boundary condi- 
tions produce inhomogeneous density distributions even 
in the absence of external fields. 

We have used the LADA approximation ((85|) . with 
(n) approximated by the parametrization for e^**' (n) 
given by (PT|) - ([M|) . to obtain the effective magnetic fields 
{/if, . . . , entering the Kohn-Sham Hamiltonian (|30p . 
Through a Jordan- Wigner transformation^^ this Hamil- 
tonian is easily diagonalized for lattice sizes up to L « 
3000. The density of up spins in the ground state is 
then calculated by diagonalizing (f5Dl) with (|95p . In order 
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to compare our LADA predictions with exact results for 
small chains, we present in table |TT] some of our spectral 
calculations at anisotropy A = —1/2. The table displays 
the lowest energies E°^°^{L, 5 + x) ^^'^ mass gaps 
(see ([68)) ) in the sectors with N = -j + up spins 
[v ~ 0,1,2) and lattice sizes L = 4 to 24. The exact 
results were obtained by a direct diagonalization of ([25| , 
and 7=1 was used in the LADA functional (|55|) . We ob- 
serve good agreement of the LADA predictions with the 
exact data, becoming better as the lattice size increases. 

As discussed in Sec. IIVI the mass-gap amplitudes of 
the finite-size corrections give the surface critical expo- 
nents x'^^ of the models (see Eq. ([681)). We verified, 
for arbitrary values of A, that the LADA functional ([551) 
with 7 = 1 and the parametrization ([9ip gives quite good 
predictions for the surface exponents. In tables IIIII and 
IIVI we show some of our estimates obtained for lattice 
sizes up to L = 512. The finite-size estimates for the 
exponent x''^ with v = \ and the ratio x"^ /x"^ are pre- 
sented in the last columns of these tables. Exact results 
in the thermodynamical limit L — > 00 are shown in the 
last line of the tables. These results were derived from 
the Bethe ansatz equations (see appendix) and conformal 
invariance (see ([5^ and (|24p ). 

It is interesting, at this point, to discuss the effect 
of the parameter 7 in our LADA functional ([85]) (or 
(|84p ). Strictly speaking, 7 is a phenomenological param- 
eter whose value was argued to be 7 « 1 in Sec. IV Al 
The results presented in table IIIII and IIVI were obtained 
by setting 7 = 1. We have also calculated the quantities 
presented in table HIT] and IIVI for other values of 7. For 
A < we got good results even from 7 = 0, where the 
nonlocal correction vanishes and the LADA functional 
reduces to the LDA. However for A > (ferromagnetic 
regime) we get non- vanishing gaps for 7 = 0. This means 
that in this regime nonlocality is important for inhomo- 
geneous density profiles and cannot be neglected. Good 
results are obtained for all values of —1 < A < 1, as long 
as 0.98 < 7 < 1.03, which is consistent with our a priori 
expectation 7 « 1. 

Let us now consider the density profiles {ni ,...,71,^,} 
of up spins on finite lattices. The Bethe ansatz for the 
finite chains, which would give the exact results for the 
density profiles, is impractical to solve. We have there- 
fore considered lattice sizes we are able to diagonalize 
numerically. In Fig we show, for L = 24 sites, an ex- 
ample of the exact and LADA prediction for the density 
profile {ni, . . . , 7124} of the XXZ chain with open bound- 
aries. The densities shown in the figure are obtained 
from the lowest eigenenergy in the sector with n = 11 
up spins and anisotropy A = — i. We see from this fig- 
ure that although the LADA amplitudes are smaller than 
the exact ones, they exhibit the same spacial oscillations. 
These spacial oscillations should decrease in amplitude 
as the lattice size increases. In Fig [6] we show the den- 
sity profiles predicted by our LADA functional ([85]) with 
the parametrization function (|9ip and 7 = 1, for sev- 
eral lattice sizes. The density profiles correspond to the 



1 spins at 



lowest eigenstate in the sector with N = ^ 
anisotropy A = — i. We have scaled the vertical and hor- 
izontal axis of Fig. [51 in order to display simultaneously 
the density profiles of distinct lattice sizes. This figure 
shows that the amplitude of the oscillations decreases as 
0(L^^), as the lattice size increases. 



0.6 



0.55 



C OAi 



0.35 



LADA 
exact 




20 



FIG. 5: (Color online) Exact and LADA predictions for the 
density profile n(i) of up spins for the XXZ quantum chain 
with open boundaries and A = — i. The profiles correspond 
to the lowest eigenenergy in the sector with iV = 11 up spins 
of the quantum chain with L = 24 sites. 



C. Periodic boundary condition with 
inhomogeneous external magnetic fields 

A physically distinct way to produce inhomogeneous 
density profiles {m, . . . , n^} of up spins is the presence 
of external magnetic fields {ft-i"*, . . . , /i^f*}- In this case, 
as in the case of open boundaries, we have to solve the 
Kohn-Sham Hamiltonian ([30|) - ([3T|) selfconsistently. 

In order to compare our results to the exact ones, we 
show in Fig. [3 the density profile obtained by using the 
LADA functional ([84|) with e^(ri) replaced by e^"(7T,) 
given by (|9ip . The density profiles of Fig. [3 correspond 
to the lowest energy in the sector with N ^ 12 up spins of 
the periodic XXZ chain with L = 2A sites and anisotropy 
A = — i. The external local magnetic fields are chosen 
to be /i™' = Si^ii- Similar to Fig. [51 the LADA densities 
have smaller amplitudes than the exact ones, but exhibit 
the same spacial oscillations. 

One of the main purposes of this paper is to present 
density functionals that are able to predict the critical 
exponents of critical chains from finite-size corrections 
to the mass gap. In order to test the LADA functional 
([84[) with ([9T|) we study the periodic XXZ chain with 
a single impurity, represented by the external magnetic 
field ft,™* = 5i^i {i = 1,...,L). The presence of such 
impurity breaks spin reversal symmetry, but we expect. 
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TABLE II: Lowest eigenenergies in the sector with N = ^ + u (u — 0,1) up spins of the XXZ chain with open boundaries with 
A = — i and L = 4 — 24. The 2nd and 4th (3rd and 5th) columns are the exact results (LADA results). The exact and LADA 
results for the gaps are shown in the 6th and 7th columns, respectively. 
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-0.733365 


-0.726786 


-0.720956 


-0.713716 
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TABLE III: Lowest eigenenergies in the sector with N = ^ + u {u — 0,1,2) up spins of the XXZ chain with open boundaries 
and anisotropy A = — i. The results were obtained by taking 7 = 1 and using the parametrization (|9Hl - (|92p in the LADA (|85|) . 
In the 5th and 6th columns we show the finite-size estimates for the surface critical exponents x"^^ and the ratios x'^'^/x"^^, 
respectively. The last line gives the expected exact results in the bulk limit {L —> 00). 
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TABLE IV: Same as table IIIll but for anisotropy value A — 
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FIG. 6: (Color online) LADA predictions for the density pro- 
files n{i). The profiles correspond to the lowest eigenenergy 
in the sector with N = ^ — 1 up spins of the XXZ quantum 
chain with L sites (L — 64,128,256), open boundaries and 
anisotropy A = — ^ . 




FIG. 7: (Color online) Exact and LADA results for the den- 
sity profiles corresponding to the ground state of the XXZ 
quantum chain with periodic boundaries and inhomogeneous 
magnetic fields. The profiles correspond to the quantum chain 
with L = 24 sites, A'' = 12 up spins, anisotropy A = — i, and 
inhomogeneous magnetic fields = Si^n (i = 1, . . . , 24). 

as i — > cxD, the leading behavior for the average mass gap 

~2Er{L,^)] (96) 

TT ^ ^ 1 

= -^WAX^" -f-o(-), (97) 

where va is the sound velocity (|45| and a;^"^ the surface 
exponents {v = 1,2,...). In Eq. ^ E^^{L,n) is the 



TABLE V: LADA predictions for the ground state energy and 
finite-size estimates of the average mass gaps (|96p for the XXZ 
chain with periodic boundaries and inhomogeneous magnetic 
fields. The results were derived for the anisotropy value A = 
— \ and inhomogeneous fields hf^"^ = 5i,L/2 (i = 1, . . . ,L). In 
the last line it is shown the expected results in the bulk limit 
(L -> oo). 
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lowest energy in the sector with N — nL up spins of 
the periodic chain with anisotropy A and L sites. In 
tablelVlwe present the finite-size estimates for the critical 
exponents a;^'' of the quantum chain with A = — i and 
lattice sizes up to L = 1024. In the last line we show 
the expected exact results in the bulk limit L — > oo. 
Clearly, the LADA functional (|84p with parametrization 
(pij) gives reliable results for the critical exponents. 

VI. CONCLUSIONS 

In this work we explored the possibility to obtain reli- 
able results for critical exponents, conformal anomalies, 
and related quantities of quantum chains from density 
functionals. Due to conformal invariance, the critical ex- 
ponents are obtained from the 0(1/ L) finite-size correc- 
tions of the mass gaps of the quantum chains in finite 
geometries. These finite-size corrections can be obtained 
approximately from DFT. In order to test our general ap- 
proach, we performed an extensive study of various LDA 
and beyond LDA functionals for the exactly integrable 
XXZ quantum chain in the critical regime — 1 < A < 1. 

Our functionals are obtained from a formal gradient 
expansion of the unknown exact functional around ho- 
mogeneous distributions of the infinite system. The first 
term of this expansion gives LDA functionals, appropri- 
ate for cases where the density distribution in the chain is 
weakly inhomogeneous, or fully homogeneous. The latter 
is the case when the quantum chain is defined on trans- 
lationally invariant lattices, with periodic and twisted 
boundary conditions (see Eqs. ([55)1 and ([55)1 ). As we 
saw in Sec. IIV[ even a simple LDA can give exact crit- 
ical exponents, but a direct evaluation of the conformal 
anomaly from these functionals gives wrong results be- 
cause the sound velocity, a nonuniversal constant, has 
distinct values for the interacting Hamiltonian and the 
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non-interacting auxiliary Kohn-Sham Hamiltonian. A 
correct LDA prediction of the conformal anomaly c = 1 , 
for any value of the anisotropy A, can be obtained only if 
we rescale the interacting and non interacting quantum 
chain by its appropriate sound velocity. This is certainly 
a problem for general quantum chains whose sound ve- 
locity is unknown. A possible solution is the use of a size- 
dependent LDA, which embodies nonlocal corrections in 
a simple way. 

In the case of open boundaries, only eigenstates be- 
longing to the sector with N = L/2 up spins have a 
homogeneous distribution on the quantum chain with 
finite size L. In this case, the surface energy of the 
quantum chain is predicted incorrectly, for the same rea- 
son discussed above for the sound velocity and confor- 
mal anomaly. However, as we showed in Sec. IIVI (see 
Eq. ((77)l ). application to a properly rescaled Hamiltonian 
gives correct results, as does use of a finite-size LDA. 

If the density distribution is inhomogeneous, our re- 
sults indicate that we must also consider at least the 
second term of the expansion (|15 p - (fT5|) of the exact func- 
tional. The LADA functional derived from this term for 
quantum chains with periodic and open boundary con- 
ditions is given by ((84)) and ([85)l . respectively. These 
functionals can be used, in principle, for a general quan- 
tum chain, by replacing e^(n) by an approximate ex- 
pression e^°™(rt) for the ground state energy per site of 
the infinite system with density n. In this general case, 
the quality of the resulting prediction of critical expo- 
nents, via evaluation of gaps, depends on the quality of 
approximations for e*^°™(n). For an exactly integrable 
chain, like the XXZ quantum chain, the energy per site 
e^(n) can be obtained exactly, e.g., by solving the inte- 
gral equation derived from the Bethe ansatz (see (jA.l[) - 
(|A.3[) in the appendix). Instead of using such numerical 
results for e^(n) we proposed, in Sec. |Vl an approximate 
parametrization e^^{n), Eq. (|9ip . containing all essential 
ingredients to furnish good estimates for the critical ex- 
ponents of the XXZ quantum chain in the presence of 
small inhomogeneities. 

Both the finite-size LDA and the LADA functional are 
nonlocal density functionals, the former depending on 
the local density and the size of the system, the latter 
depending on the local and the average density. Unlike 
the usual (infinite-size) LDA, these functionals also de- 
pend, in a known way, on the boundary conditions (open, 
periodic, or twisted). As a consequence of this nonlocal- 
ity, the finite-size LDA and the LADA functional each 
cure certain defects of the LDA. For finite homogeneous 
chains, the LDA reproduces the correct critical expo- 
nents, but predicts a wrong conformal anomaly, whereas 
the finite-size LDA correctly reproduces both. For in- 
homogeneous chains, the LADA functional yields better 
energies and gaps than the LDA. For finite homogeneous 
chains, the LADA functional reduces to the LDA. The 



combination of both improvements, i.e., the construction 
of a finite-size LADA is a promissing project for the fu- 
ture. 

Another interesting question concerns possible exten- 
sions of the LADA functional (j84p - ([85l) to other quantum 
chains in the c = 1 universality class. Consider a general 
homogeneous model on this universality class. The in- 
troduction of a homogeneous magnetic field (or chemical 
potential) h{n) that couples to the density operator fii 
will produce a ground state with energy per site e(n) and 
homogeneous density n. In order to fix the density n, this 
magnetic field should satisfy h{n) = The critical 

exponents with such external field are given by Eq. 
with X = x{n), depending on the particular density fixed 
by h{n). Besides x{n) and h{n), the quantum chain is 
also characterized by the non-universal sound velocity 
v{n) that fixes the scale of the momentum-energy dis- 
persion relation. The results coming from conformal in- 
variance for this class of models give us a generalization of 

Eq. (|48| for arbitrary densities, i.e., = AT:v{n)x{n). 

In order to obtain the density distribution {n^} when 
inhomogeneous external fields or boundaries are added to 
the homogeneous system we may proceed in two ways. 
If the local density never exceeds unity, the XY model 
pop can be used as the auxiliary non-interacting (Kohn- 
Sham) system. If this is not the case, the only remain- 
ing possibility is the direct use of the approximate func- 
tional for F] in the Euler equation ([6]). The first case is 
preferable, since the LDA functional is not exact and the 
auxiliary non-interacting chain will give a more precise 
evaluation of the kinetic energy contained in . In both 
approaches, the necessary quantities for the evaluation of 
the density are h(n) and 47rw(rt)x(n). 

To conclude, let us briefly discuss some possible fu- 
ture developments of the present work. The functional 
we produce gives us mass gaps of finite chains that pro- 
duce reasonable, in some cases even exact, estimates for 
the critical exponents. A second possibility to obtain 
the critical exponents is from the power-law decay of the 
Friedel density oscillations due to the presence of im- 
purities on the quantum chaini^ Tests with the LADA 
functionals (|84p or (|55|) for the XXZ chain with open and 
periodic boundary condition indicate poor results in this 
case. This means that the correct decay of the Friedel 
density oscillations can only be obtained by considering 
higher order terms in the expansion p5p - (fT5|) . This is 
certainly an interesting point for the future. 
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APPENDIX: INTEGRAL EQUATIONS DERIVED 
FROM THE BETHE ANSATZ SOLUTION OF 
THE XXZ QUANTUM CHAIN 

The ground state energy per site, e^(rt), of the XXZ 
chain, in the bulk limit, can be obtained from the Bethe 
ansatz solution of the model (see e.g. Ref. [SC 

For — 1 < A < 1 , by setting A = — cos 7 



cos 7 



sin^ 7 



cosh rj ~ cos 7 



R{T^)dr^, (A.l) 



where the parameter A fix the density n through 
R{ri)dr] — n. 



(A.2) 



The function R{ri) is obtained by solving the integral 
equation 
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cosh(?7 — rj') — cos 27 



cosh 77 — cos 7 
R{rf)dri' 
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